Violating conformal invariance: Two-dimensional clusters grafted to wedges, cones, 

and branch points of Riemann surfaces 
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Lattice animals are one of the few critical models in statistical mechanics violating conformal 
invariance. We present here simulations of 2-d site animals on square and triangular lattices in 
non-trivial geometries. The simulations are done with the newly developed PERM algorithm which 
gives very precise estimates of the partition sum, yielding precise values for the entropic exponent 
9 (Zjq ~ fj,^ N~^). In particular, we studied animals grafted to the tips of wedges with a wide 
range of angles a, to the tips of cones (wedges with the sides glued together), and to branching 
points of Riemann surfaces. The latter can either have k sheets and no boundary, generalizing in 
this way cones to angles a > 360 degrees, or can have boundaries, generalizing wedges. We find 
conformal invariance behavior, ~ l/a, only for small angles {a <^ 27r), while 6 ~ const — a/2TT for 
Q S> 2n. These scalings hold both for wedges and cones. A heuristic (non-conformal) argument for 
the behavior at large a is given, and comparison is made with critical percolation. 



Lattice animals (or polyominoes, as they are sometimes 
called in mathematics 0) are just clusters of connected 
sites on a regular lattice. Such clusters play an important 
role in many models of statistical physics, and are consid- 
ered as the standard model for randomly branched poly- 
mers The difference between the animal model and 
other cluster models such as percolation is that the clus- 
ters appear with non-trivial weights in the latter, while 
every cluster with the same number N of sites has the 
same weight in the animal ensemble (this gives, more 
precisely, site animals; in this letter we shall not con- 
sider other models in the animal universality class, such 
as bond animals or lattice trees Q). Moreover, lattice 
animals are one of the few critical models in statistical 
mechanics violating conformal invariance 0, 1^- How- 
ever, the deeper reasons for this violation as well as the 
consequences are still poorly understood. 

On the one hand, the animal problem is similar to other 
models of statistical physics in allowing a field theoretic 
formulation ^ and in showing anomalous scaling laws in 
space dimensions d < dc = 8, where dc is called the upper 
critical dimension. In particular, the number of animals 
with precisely N sites attached to a given fixed site (the 
"partition sum" ; notice that we count here shifted ani- 
mals as different) scales for large N as 



(1) 



where fj, is the non-universal growth constant and 9 is 
a critical exponent which is independent of the lattice 
type, but depends on the global geometry of the lattice. 
Similarly, the gyration radius scales as i?Ar ~ N'^ where 
the Flory exponent v is universal and also independent 
of the geometry of the lattice. 

On the other hand, the statistics of lattice animals 
has a number of rather unusual features. In particu- 
lar, as mentioned above, it is not conformally invariant 

. Conformal invariance gives rather strong constraints 
in two dimensions. First of all, it gives restrictions on 
critical exponents. The fact that the Flory exponent v 
is not exactly known for 2-d lattice animals (while it is 



known exactly for unbranched polymers and most other 
2-d models), is a consequence of the lack of conformal 
invariance. 

Secondly, and more closely related to the present work, 
conformal invariance gives strong constraints on the en- 
tropic exponent B in non-trivial geometries. The most 
thoroughly studied of such geometries are wedges and 
cones. A wedge is a part of the plane bounded by two 
straight lines which intersect at an angle a. A cone is 
basically a wedge where the two boundaries are glued to- 
gether. A wedge can be mapped onto a half plane by a 
conformal map, while a cone can be mapped onto a punc- 
tuated plane. This implies that for conformally invariant 
theories the entropic exponent is linear in 1/a 



(2) 



For models in which one can study single clusters (such 
as self avoiding walks, percolation, or lattice animals) 
this applies to clusters grafted at the tip of the wedge 
or cone, respectively. While this equation was checked 
for self avoiding walks (linear polymers) 0,0], it was 
indeed found not to hold for lattice animals |a, Q • But 
the numerical results given in the latter papers were not 
sufficient to suggest any alternative behaviour. 

In the present letter we apply the PERM (Pruned- 
Enriched Rosenbluth Method) strategy which was re- 
cently adapted to lattice animals It is a recursively 
implemented sequential Monte Carlo method with re- 
sampling which starts off by growing percolation clus- 
ters, re- weighs them according to the animal ensemble, 
and applies cloning and pruning to achieve approximate 
importance sampling. It is the most efficient algorithm 
for simulating animals and lattice trees known today. In 
particular, it provides very precise estimates of the par- 
tition sum, which then allows to estimate Q by means of 
Eq.Q. In the following, the maximal sizes of animals 
varied between iVmax — 1000 and iVmax = 4000. 

The results presented below are based on simulations 
on the square and on the triangular lattices. While the 
growth constant for the square lattice was taken from 
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TABLE I: Entropic critical exponents for 2-d lattice animals 




FIG. 1: Typical cluster for a wedge with a — tt/3 on a trian- 
gular lattice; N = 3500. 



0, we had to perform new simulations on the full lat- 
tice to obtain it for the triangular lattice. Wc obtained 
In/itriang = 1.6454139(18) (using the fact that 6* = 1 on 
a full 2-d lattice). 

Wedges on the square lattice are most easily obtained 
by placing one edge along one of the coordinate axes (say 
the positive x-axis), and taking the angle a such that 
tana = n/m with n and m being integers. Alternatively, 
one can place the x-axis along the center of the wedge, 
and use tan(a/2) = n/m ("symmetrical wedge"). While 
these constructions cannot be used e.g. for wedges of 30, 
60, or 120 degrees, the latter can be obtained for trian- 
gular lattices, see Fig. ^ We checked that the standard 
and the symmetrical wedge gave the same 9 for 90 degrees 
(within error bars), and that the square and triangular 
lattices gave the same results for 180 degrees. Wedges 
with 360 degrees are obtained by excluding a single half 
line. 

As we said, cones are obtained by gluing together the 
two edges of a wedge. This can of course only be done if 
the lattices agree along the two edges and if this does not 
introduce a line of defects, which strongly restricts the 
possible angles. For the square lattice, only a = 90, 180, 
and 270 degrees are possible, while 60, 120, 240, and 
300 degrees are also possible for the triangular lattice. 
Enumeration data for cones with other angles are given 
in j^, but it is not clear how they were obtained. In any 
case we checked that defect lines per se have an effect 
on 9, by simulating animals grafted to points on a defect 
line. 

To simulate angles larger than 2tt, we used multi- 
sheeted Riemann surfaces with branch points at the ori- 
gin where the animal is grafted. A branch point where k 
sheets meet is essentially a cone with angle 2fc7r. Wedges 
with a > 27r are then obtained by cutting out the corre- 
sponding domain from the surface. 

There are a number of implementation details which 
have to be specified. For instance, one has to specify for 
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wedges whether a site on the boundary can be occupied 
or not. For a cone one might specify whether the two 
boundaries are identified, or whether they are one lat- 
tice unit apart. Finally, for Riemann surfaces the site at 
the origin can have M neighbours {M is the coordination 
number. A/" = 4 for the square lattice and A/" = 6 for 
the triangular one), all of them on one sheet. Or it can 
have kN neighbours, occupying all k sheets. We checked 
in each case several of these alternatives. As expected, 
they gave different results for finite N , but they led to 
the same scaling behaviours. 

Our final results are given in Tables 1 and 2, and are 
also shown in Fig. |2] Notice that there are two values 
which are exact: 9 = \ for cones with angle 27r, and 9 = 2 
for wedges with a = tt j^] . For a ~ tt our data agree very 
well with those of but for smaller a the latter data 
seem to be systematically too low: The value 9 = 5.5±0.1 
for a = 30 degrees cited in e.g., is seven standard de- 
viations below our value 6.204 ± 0.008. A detailed com- 
parison with the cone data of is less straightforward. 
For those angles where cones without defect lines exist, 
the agreement is excellent. For those where there should 
have been defect lines (127°, 143°, 233°), the data of [3 
seem to be too high by two to four standard deviations 
(when compared to smooth interpolations of our data), 
as should be expected from the above discussion. 

Both data sets (wedges and cones) confirm that Eq. Q 
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TABLE II: Entropic critical exponents for 2-d lattice animals 
grafted to the tip of a cone with angle a. 
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FIG. 3; Values of the exponent r — 2 for critical percolation 
clusters, grafted on wedges with angle a. The points (whose 
error bars are smaller than the point size) are from simulations 
for site percolation on the square lattice, the straight line is 
the theoretical prediction r — 2 = 167r/9lQ. 
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FIG. 2: Values of the entropic exponent 9 for wedges (upper 
curve) and cones (lower curve) with angle a, plotted against 
a. Values for a > 2tt were obtained using multi-sheeted Rie- 
mann surfaces. Statistical and systematic errors are much 
smaller than the sizes of the symbols. 



does not hold. But the wedge data (Table 1 and Fig. ^ 
show very clearly that Eq. (|2Jl does hold asymptotically for 

9{a) w 1.35 + 6wodgo/a with 6wodgc = 2.543 ± 0.020. 

(3) 

A similar scaling is also compatible with the cone data, 
although there the error on b is much larger since we could 
not go to sufficiently small angles: 6conc = 1-4 ± 0.1. On 
the other hand, both data sets indicate that 6 increases 



linearly for large angles, 

9{a) a for a — > oo. (4) 

Moreover, the coefficient of proportionality seems to be 
exactly the same in both cases, 

lim --^ = (1.00 ±0.03). (5) 

Q— ►oo a ZTT 

In order to understand this behaviour heuristically, we 
now compare lattice animals (which are in the univer- 
sality class of subcritical percolation) to critical perco- 
lation. If we use a standard cluster growth algorithm 
like the Leath algorithm 01 or the depth-first algorithm 
of Swendsen and Wang [T^l to grow critical or slightly 
subcritical percolation clusters, the probability to reach 
a size > N decreases like a power 111 , 

PN-^N^-^f{{p,-p)N^). (6) 

Here, p is the wetting probability, pc its critical value, 
T and a are critical exponents, and f{x) decreases ex- 
ponentially for a; — > oo. The critical exponents can be 
related to other, more standard, exponents, e.g. |llj | 

r = 3 — 7cr. (7) 

The ansatz Eq.ljHJ holds both for clusters grown in the 
bulk, and for clusters grown near a surface. In the latter 
case, Pc and the exponent a = 36/91 are the same as in 
the bulk (and independent of the shape of the surface) , 
while T and fix) do depend on the surface. In particular 
it is known jj| that 7 = 25/12 for 2-d clusters attached 
to a plane wall, giving r(a = tt) — 2 = 16/91. 

From conformal invariance we expect that r(a) is a 
linear function of 1/a. On the other hand, we expect 
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that r — * 2 for a — > cxD. The reason is very simple: 
For a ^ oo, the chances that the growth will stop at any 
finite N will go to zero, since there are ever more possible 
directions for growth. Thus we expect 



From Fig. 13 we see that this is in excellent agreement 
with simulations of clusters starting at the tip of a wedge. 
Notice that the angles shown in Fig. |31 extend up to Stt, 
verifying thereby that branch points of Riemann surfaces 
can be treated like tips of wedges. 

For lattice animals one could also try to use Eq.®, 
this time with p pc, but this would not lead to any 
useful prediction. Thus one has to proceed differently. 
Our main assumption is that clusters will grow essentially 
into an angular region of size Aa = 0(1). Much larger 
angular ranges will also occur, but only with very low 
probability. For a :s> Aa, i.e. k = a/Aa :s> 1, one has 
thus essentially k independent clusters. If one assumes 
that all these subclusters have roughly the same size, one 
expects 

Zr,ika) « [Z^/kiat. (9) 

Together with Eq. this gives 6* oc a for large a. This 
argument explains also why the proportionality constant 
is the same for wedges and cones, but it does not explain 



its numerical value of 1/2tt, i.e. it does not explain why 
6 decreases by exactly one unit when a increases by 27r. 

If lattice animals are not conformally invariant but in 

some way "covariant" , one might expect a simple analyt- 
ical formula for 9{a). We therefore tried to find such fits. 
Simple ansatzes like 6 = a + ba + c/a were not successful. 
The simplest acceptable fits were obtained with Fade ap- 
proximates of the form 9 — {a+ba+ca^ + da^) / {a+ea^). 
But the coefficients a to e looked rather uninspiring, and 
such an ansatz seems already too complicated for being 
"natural" . 



In summary, we have presented very high statistics 
simulations of clusters grafted to the tips of cones and 
wedges. For critical percolation clusters these simula- 
tions were in full agreement with predictions from con- 
formal invariancc. But for lattice animals (subcritical 
percolation clusters) they agreed with the conformal in- 
variancc behavior only in the limit of small angles. For 
large angles another simple behaviour was found and ex- 
plained by assuming distant angular regions to be es- 
sentially independent. Our results clearly show conse- 
quences of the violation of conformal invariance in the 
lattice animals model. Nevertheless, the results also sug- 
gest that for small angles conformal invariance may still 
hold in some approximate way or that some generalized 
invariance might exist. 
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